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Abstract 

In order to investigate the dispersion of the spectral break energy in gamma- 
ray bursts, we simulate internal shocks, including effects of shell-splitting after shell 
collisions and the Thomson optical depth due to electron-positron pairs produced 
by synchrotron photons. We produce pseudo observational data and estimate the 
break energy. If the distribution of initial Lorentz factors of shells has only one peak, 
many pulses with a break energy much lower than the typical observed value should 
be detected, while the effect of the Thomson optical depth reduces pulses with a 
break energy of > 1 MeV. If the Lorentz factor distribution has multiple peaks, the 
distribution of the break energy can be consistent with observations. 

Key words: gamma rays: bursts — radiation mechanisms: non-thermal — shock 
waves 



1. Introduction 



Gamma-ray bursts (GRBs) are short 10 s) bursts of low-energy 7-rays. One of the 
most important characteristics of the GRB spectra is the existence of a typical break energy 
scale. The observed spectra of GRBs are approximated by a broken power law and the photon 
number spectra are approximated as ~ above the break energy and ~ below that. The 
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standard scenario for producing GRBs is the dissipation of the kinetic energy of a relativistic 
flow by relativistic shocks (see, e.g., a review by Piran 1999). The rapid time variabilities 
observed require that the GRB itself must arise from internal shocks within the flow. The 
radiation is emitted by relativistic electrons in the shocked region. The observations are well- 
described by synchrotron emission (Cohen et al. 1997; Wijers, Galama 1999). Kobayashi et 
al. (1997) shows that the internal shock can reproduce the temporal structure of GRBs. A large 
fraction of the kinetic energy could be converted to radiation by the internal shocks (Kobayashi, 
Sari 2001). 

The internal shock model is roughly represented by colliding shells. The collisions are 
considered to occur within a wide range of radius r. The energy efficiency argument requires 
a large dispersion of the initial Lorentz factor F (Kobayashi, Sari 2001; Beloborodov 2000). In 
spite of the large dispersions of r and F in the model, the apparent clustering of the break 
energy of the GRB spectra in the 50 keV-1 MeV range is reported in the BATSE observation 
(Preece et al. 2000). It is strange that pulses in a burst tend to have similar break energies. 
The break energy distribution is in agreement with a log-normal distribution. This distribution 
seems to suggest small dispersions of r or F. 

In this paper, based on numerical simulations, we examine whether the standard internal 
shock model can reproduce the narrow distribution of the break energy. Guetta et al. (2001) 
showed that the internal shock model can reproduce the typical break energy. However, the 
dispersion of the break energy has not been explained so far. In our simulation, as Guetta et 
al. (2001) did, we include the effect of the Thomson optical depth due to pairs produced 
by synchrotron photons. We treat the pair optical depth more precisely. In addition, our 
improvements in the numerical simulations over those in Guetta et al. (2001) take into account 
shell splitting and the spectral energy band in the observation. From our simulation we obtain 
some restrictions on GRB models. In section 2 we explain our simulation method. In section 
3 the numerical results are given. In section 4 we summarize our results. 

2. Method of Simulation 

In order to examine the distribution of the break energy, we consider a spherical wind 
consisting of shells. Each shell is characterized by the Lorentz factor Fj, mass Mj, width 
Wi, and radius from the center r^. Following the evolution of the shells, there will be numerous 
collisions between different shells. For each collision we calculate the break energy, optical 
depth, and flux. In order to simplify our model we idealize the situation. Some approximations 
are adopted in our formulation, as discussed below. 

2.1. Break Energy 

In this subsection we explain our method to obtain the break energy. Let us consider 
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shells with different Lorentz factors, x = Fr/Fs > 1. The relative Lorentz factor is Frci ~ 
(x + l/x)/2. The rapid and slower shells are denoted by subscripts r and s, respectively. 
We assume that the widths of the shells are comparable in the interstellar medium (ISM) rest 
frame as Wi = W. Under these assumptions, we can obtain the ratio of the baryon number 
density, Uj-fus. For equal mass shells {Mi = const.) rij-fus = l/x, while rij-fng = for equal 
energy shells (MjFj = const.). In this paper, n and internal energy density e are measured in 
the fluids' rest frame. 

When the rapid shell catches up with the slower one at some radius r, the forward and 
reverse shocks form. The shock conditions and the equality of pressures along the contact 
discontinuity (Sari, Piran 1995) yield 

(Ff-1)(4Ff + 3) ^ 

(Fr-1)(4Fr + 3) ns' 
where Fp and Fr are the Lorentz factors of the relative motion between the regions separated 
by the forward shock and by the reverse shock, respectively. The equality of the velocities 



along the contact discontinuity gives Fr = FpFrei — y (F| — l)(F^pj — 1). In the case of equal 
mass shells or equal energy shells, the solution for Fr depends on only the ratio of the Lorentz 
factors X. The internal energy density in the shocked regions is given by 
(Fr-1)(4Fr + 3)M,c^ 

where Mr is the mass of the rapid shell. 

We now consider the synchrotron emission from shocked shells. Since a rapid shell has 
a larger energy and a lower number density in the cases of equal mass and equal energy, the 
average energy per one electron in the reverse shock is larger than that in the forward shock. 
The reverse shock emission is more luminous and harder than the forward shock emission. 
Therefore, we consider only the emission from the reverse shock. The shock is assumed to 
accelerate electrons in the shell material into a power-law distribution, A^(7e) oc 7^^(7e > 7c,min)- 
Assuming that constant fractions, and ee, of the internal energy go into the magnetic field 
and the electrons, respectively, one finds that the magnetic field and the typical random Lorentz 
factor of electrons are given by B"^ oc e^e and 7e,min oc ec(FR — 1), respectively. 

The synchrotron process is a very efficient radiation process. With the strong magnetic 
field required to produce the observed gamma-ray, the synchrotron cooling time is very short 
compared to the dynamical time of the shock. In this fast cooling case (Sari et al. 1998), the 
photon number density distribution at the fluid rest frame is given by 

dn^{e) _ p-2 ecC j {e/e^Y'^''^ for 5 < ^p, 

de ■ j (,/,p)-(^+2)/2 for5>5p, 

where the break energy is the typical energy of synchrotron photons emitted by electrons of 



7e,min- Siucc the shocked fluid is moving with the Lorentz factor, Fm ~ Fi.(Fr — — 1), the 
observed break photon energy is 
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£°^^ = 610e2F(x)r,,2v/s7A keV, (4) 

where F{x) = x^/^ (Pj^ -1)5/2 ^4^^ ^ 1/2 _ _ 1 ) . It behaves as ~xV4 at x > 1. We have 
scaled the parameters as €3-1 = 6^/0.1, M48 = ergs, ri3 = cm, W7 = W/IO'^ 

cm, and rs^2 = Ts/lOO. We adopt p = 2.5 here and hereafter. Since some parameters appear 
in the following formulae in the same combinations, we have defined two variables, a surface 
density S = M^s/rf^ and an effective comoving shell width A = W7rs,2/^B,-i for convenience. 
In our simulation we neglect the effect of the cosmological redshift on the break energy. 

2.2. Optical Depth Due to Pairs 



Shells are initially optically thick due to Thomson scattering by electrons associated with 
baryons in the shell. For a large optical depth r, the radiative cooling time can be estimated 
as ~ tI/c, where / = WTj-f^AT^ + 3) is the shell width in the comoving frame. Comparing 
this with the cooling time due to the shell spreading ~ //c, the radiated energy is negligible 
for r ^ 1. Then, all of the internal energy is transformed back to kinetic energy (Kobayashi, 
Sari 2001). Thus, we could not observe the emission from it, until it comes from outside of the 
photosphere, where the optical depth tt = 2crTM/47rr^mp ~ 0.7S becomes unity. This is a well- 
know result, but the estimate may undergo a significant change in the internal shock model if 
we take into account pairs produced by the synchrotron photons (Guetta et al. 2001). In the 
fluid rest frame, the break photon energy Sp is much smaller than the electron mass. However, 
the photon distribution extends to high-energy as a power law. The pairs caused by the high 
energy photons may contribute significantly to the optical depth. 

When a photon of energy e interacts with a photon of e' at an incident angle 6, the cross 
section of the pair creation is written as a± = crTf{y) (Berestetskii et al. 1982), where 



f{y) = ^{i-y') 



(3-y^)ln\±y-2yi2-y') 



(5) 




The dimensionless value y is defined by = 1 — {2mlc'^)/[ee'{l — cos6')]. The optical depth to 
pair-creation is given by 

T^y{e) = f {l-cose)dQ I ^^2,4 ^g' Tj^n ^77(g'^g)^ (6) 

(p-l)(p + 4)£ 1^^'"^' 
where 

C{p)^ f'dy{i-yY'~'yfiy)- (8) 

Jo 

The value of C{p) is not sensitive to p, and C(2.5) ~ 0.075. Since we use the power-law 
distribution dn^/de the integration, the optical depth is overestimated for photons 

with high energy e ^ 2m^c"^/£:p, which interact with photons mainly in the low-energy portion 
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dn^/de oc e"^/^. However, the exact value is not important for the high-energy photons when 
we estimate the number of pairs below, because the corrected optical depth is also large 
enough to annihilate all of them. For p = 2.5 we obtain 

r,,{e) = 3.3 el^'Gix)!:^/'' A^'/^ [e/m,c^f\ (9) 

where G{x) = (Fr - l)i3/8(4rR + 3)^/Vxi/^ Although the optical depth to pair-creation for 
photons with the typical energy of Ep <^ meC^ is very small, it does not necessarily imply that 
the number of pairs is negligible for Thomson scattering. Considering that photons with 
energy e interact mainly with those with energy ~ 2m^c^/e to produce pairs, we can estimate 
the optical depth due to the pairs as 

r± ~ 2aT/ r de^^^il - e-^-). (10) 

JV2m^c2 de 

If we consider only the effect by the electrons accompanied by baryons, the expanding 
ejecta becomes optically thin when the surface density S oc decreases down to ~ 1. However, 
the pairs produced by the synchrotron photons generally increase the optical depth signifi- 
cantly. Then, the ejecta is required to expand beyond S ~ 1 to radiate. In figure 1 we plot the 
value of the surface density Em, at which the ejecta becomes optically thin, as a function of 
X. In order to plot figure 1, we have numerically solved equation (1) and integrated equation 
(10). If we describe the system as an inelastic collision between two masses, and half of the 
internal energy is converted to pairs as Guetta et al. (2001) assumed, Em is estimated to be 
~ O.OlTe'^-y/a*/ (1 + x — 1\fx). This approximation generally gives a smaller value. 

Since the break energy e°^^ is proportional to S^/^, a collision occurring with larger E 
radiates harder photons. Then, the hardest synchrotron photons are emitted just above the 
radius where S = Em for given values of the parameters x. A, and Fg. We define the maximum 
break energy Eyi as the break energy emitted from the radius where S = Em- The value of £'m 
is obtained from equation (4) with S = Em- Figure 2 depicts -Em as a function of x. We can see 
that the photon pair-creation effect makes Ey^ significantly smaller (solid lines) compared to 
the case without the effect (dotted lines). According to figure 2, it is difficult to obtain a break 
energy of more than 1 MeV for x < 20. The upper limit of the break energy is determined by 
the Thomson optical depth due to pairs, and it is close to the typical observed break energy. 
This result agrees with Guetta et al. (2001). 

In the parameter region of our interest, the critical photon energy Sr, at which the optical 
depth r^^ = 1, is larger than ^/2mec'^. Then, we obtain the approximation, 

r± ~ 300elG\x)J:^/^A-^/\ (11) 

where we neglect a logarithmic factor of e^. We adopt this approximation in our simulation. 
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Fig. 1. Plot of Em (solid line) against x. The thick and thin lines are for equal-mass and equal-energy 
cases, respectively. The dotted lines represent the values without pair-creation. For equal mass the value 
is constant as Sm ~ 1-4. Here, we assume A = 1 and = 0.5. 



10000 E r 




Fig. 2. Plot of Em (solid line) against x. The dotted lines are values obtained without the process of 
pair-creation. The thick and thin lines are for equal-mass and equal-energy cases, respectively. Here, we 
assume A = 1, Ec = 0.5, and rs.2 = 1- 
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2.3. Sampling of Break Energy 

Setting the initial distributions of the Lorentz factor Fj and separation between shells 
Li, we follow the evolution of the shells until there are no more collisions, i.e., until the shells 
are ordered by increasing value of the Lorentz factors. For each collision, using equations (1), 
(4), and (11), we estimate the break energy and the optical depth due to electron-positron 
pairs. We approximate the radiation energy as 



as conventionally assumed (Kobayashi, Sari 2001). After a fraction of the internal energy is 
emitted, the shells will split, transforming the remaining internal energy back to kinetic energy. 
In our simulation, the process of the shell splitting is basically the same as in Kobayashi, 
Sari (2001). The introduction of shell splitting in the simulation increases the number of 
collisions and the energy efficiency. In this paper, for simplicity, we assume that each mass and 
width of the shell are conserved before and after its collision. If r-t > 1, we force ec = 0. Then 
the collisions of shells for r± > 1 are similar to the perfectly elastic collisions of pool balls. 

In order to include the effect of the spectral energy band in observations, we produce 
pseudo observational data from our simulation. Using the spectrum of equation (3), we estimate 
the flux between 20 and 2000 keV from the radiation energy and break energy. We neglect the 
effect of the spectral sensitivity in BATSE instruments. The effective area changes only by a 
factor of two in most parts of the spectral energy band, so that the effect of the sensitivity 
can be negligible. From the estimated flux we write light curves by the same method as in 
Kobayashi et al. (1997). The light curves are superpositions of many pulses emitted from all 
collisions. The light curve peaks are identifled with the peak-finding algorithm described by 
Li, Fenimore (1996). A peak time Tp is identified, if the peak photon fiux Cp, photon fiuxes Ci 
(Ti < Tp), and C2 (T2 > Tp) satisfy Cp — Ci,2 > A'var-y/C^ and there are no time bins between 
Ti and T2 with photon fiux higher than Cp. We adopt A^var = 0.3. The results do not strongly 
depend on A^^ar (Spada et al. 2000). Each peak is identified as an "observable pulse" in our 
simulation. The light curve valleys are identified as the minima between two consecutive peaks. 
We divide the light curves in time into some regions by valleys and identify them as the observed 
duration time of each pulse. One observable pulse may be the composition of multiple pulses 
emitted from different collisions. For each temporal region, by attaching a fluence weight, 
we average the break energies of all pulses arrived during the period. We adopt the average 
break energy as the "observable break energy" for each temporal region. This treatment is 
different from the method used in an actual observation. In our simulation, however, one pulse 
greatly overwhelms other dim pulses in most cases. Therefore, our method is harmless when 



ee[M,(rr-r:,) + M,(r,-r'j], 



(12) 



where 




(13) 
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Fig. 3. Distribution of the break energy for all collisions. The effect of the Thomson optical depth is 
neglected. 

determining the break energy. 
3. Results of Simulation 

3.1. Continuous T- Distribution 

In our simulation, the initial shells are assumed to have equal mass M and width W. The 
distributions of the initial Lorentz factor Fj of each shell and initial separation Li between two 
consecutive shells are determined by the character of the central engine on which we have little 
information. In our simulation we assume that Fj is to be distributed uniformly in logarithmic 
space between logF^ and logFM (the number distribution of shells is proportional to 1/F). The 
initial separation Lj is also assumed to distribute in the same way between logVT and logLM- 
We perform 100 simulations with = 50, F,^ = 30, Fm = 3000, W/c = 10"^ i^^^^ ^ i 
ee = 0.6, eB = 0.1, and a total initial kinetic energy of E^so = 10^^ erg. First we present the 
results, neglecting the effect of t±. Figure 3 shows a histogram of the break energy for all 
collisions, neglecting the spectral band in the observation. The break energies widely distribute 
from the optical range to 10 MeV. This result apparently contradicts the observation. We 
then take into account the Thomson optical depth. As is shown in figure 4, the break energies 
above ~ 1 MeV are suppressed, as we have speculated in subsection 2.2. However, there remain 
many pulses with a low break energy. The effect of the Thomson optical depth reduces only 
high-break pulses. 

However, we can not observe all of the pulses emitted from all collisions. A dimmer 
pulse may be overwhelmed by other brighter pulses. Since a dimmer pulse is expected to have 
a lower break energy, actual observations reduce the number of pulses with a low break energy. 
In addition, a large energy fraction of emission with a very low break energy may be out of 
the spectral energy band in the observation. Therefore, we mimic the BATSE observation by 
producing light curves from the simulations, as explained in subsection 2.3, and define the 
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Fig. 4. Distribution of the break onoro-v fnr all mllisinns I'whifp'l and nhsorvahlr nijlses (black). The effect 

of the Thomson optical depth is t 
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Fig. 5. Same as figure 4 for observable pulses. Here we include the cosmological redshift effect, assuming 
z = f .5. The solid line is the BATSE observation. There are many other pulses outside of this energy 
band. They should be counted as pulses whose spectra have no break energy. 

observable pulses and their break energies. In figure 4, we show the result obtained from the 
pseudo observational data. Samples with lower break energy diminish drastically. In figure 5 we 
magnify the result around the BATSE band. Although the peak of the distribution is around 
the lowest limit of the BATSE band, the dispersion is significantly larger than the observation. 
There are many pulses with a break energy of < 20 keV, which should be observed as "no break" 
pulses. In the spectral energy band only high-energy tails of bright spectra of such pulses are 
detected, and they make significant peaks in the light curve. 

Regarding the observable pulses in figures 6 and 7, we show the relations between the 
break energy and the collision radius r, and between the break energy and the ratio of the 
Lorentz factors x, respectively. The collision radius and ratio x are obtained averaging in the 
same way as the break energy. According to these figures, the break energy is mainly determined 
by the ratio x rather than r. Collisions with a smaller x produce lower break energies, which 
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Fig. 6. Relation between the break energy and collision radius. 
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Fig. 7. Relation between the break energy and the ratio x. 

cause a large dispersion of the break energy. The lower limit of the collision radius in figure 6 
is determined by the Thomson optical depth. 

In this model, the fraction of the emission energy to the total kinetic energy is 16.9 ±2.0%. 
As Kobayashi, Sari (2001) claimed, the introduction of shell splitting increases the energy 
efficiency compared with that of Guetta et al. (2001) (typically a few percent). In the pseudo 
BATSE band, the energy efficiency is 2.9 ± 0.9%. About four-times radiation energy in the 
BATSE band is emitted outside of the BATSE band. 

3.2. BimodalT -Distribution 

We have simulated for other types of distributions of F. However, as long as the F- 
distribution is continuous and has only one peak, the results do not change basically. The 
models predict many "no break" pulses. Let us consider a F-distribution that has a maximum 
at F = Fp. If we randomly choose two shells in this distribution, the most probable Fs of two 
shells are around Fp. Therefore, we cannot avoid numerous collisions with x ~ 1, which leads 
to low break energies. Of course, a too small x leads to too dim emission. However, emission 
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Fig. 8. Distribution of the observable break energy for the bimodal distribution. The effect of the 
Thomson optical depth is taken into account. Here, we include the cosmological redshift effect assuming 
z = 1.5. The soUd hne is the BATSE observation. 

from a marginal value of x (~ 5) is luminous enough in spite of the small break energy. 

If the F-distribution has two peaks, the probability function of x may have a maximum 
at X ^ 1. As one example, we make a simulation for a bimodal F-distribution: one half of the 
shells have F = 30 and the other half have F = 3000 initially. The parameters are common to 
those described in the former subsection. The result is shown in figure 8. Since collisions with 
F = 30 and F = 3000 are dominant events in this case, break energies distribute narrowly. The 
dispersion of the break energy is consistent with the observation. It is trivial that we can adjust 
the peak of the distribution to the observation more closely by changing the parameter eg. 

The fraction of the emission energy to the total kinetic energy is 65.1 ±6.6%. In the 
pseudo BATSE band, the energy efficiency is 18.7 ±2.5%. The bimodal distribution is advan- 
tageous to the energy efficiency too. 

4. Conclusions and Discussion 

Following the standard scenario of GRB, we simulate internal shocks, including the 
effects of shell splitting and the Thomson optical depth due to electron-positron pairs produced 
by synchrotron photons. We produce pseudo observational data and estimate the break energy. 
The effect of the Thomson optical depth reduces the pulses with a break energy of > 1 MeV, 
which is consistent with Guetta et al. (2001). The shell-splitting effect increases the energy 
efficiency. However, many "no break" pulses should be observed in the case of a one-peak, 
continuous F-distribution. Even if we alter some assumptions (shell splitting, equal shell mass, 
etc.) in our model, the qualitative result is basically the same. Within our simple method, 
"no break" pulses are unavoidable, though we take into account the spectral energy band in 
the observation. Collisions with a small ratio x cause "no break" pulses. Preece et al. (2000) 
insists that the small dispersion of the break energy is not due to the observational selection 
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effects. However, we should henceforth consider the possibihty of selection effects on the BATSE 
observation. 

The F-distribution is determined by the character of the central engine. If the initial F- 
distribution is discrete, or has multiple peaks, like the bimodal distribution, the distribution of 
the break energy can be consistent with the observation. The bimodal distribution is favorable 
for the energy efficiency. However, we do not know if such a distribution is realistic or not. 

Judging from figure 7, if we choose pulses emitted from collisions with x ^ 10, most of 
the observable break energies are in the BATSE band. Therefore, if some microscopic processes 
prohibit collisions with x ^ 10 from emitting photons, the small dispersion of the break energy 
can be reproduced. One of the candidates for such processes is in the electron acceleration 
mechanism, which is not well understood. In the standard GRB model, a large fraction of the 
kinetic energy carried by protons is efficiently converted into that of relativistic electrons in the 
shocked region. However, this premise has not been proven theoretically. It is apparent that the 
Coulomb interaction cannot transport the internal energy of heated protons into electrons to 
achieve energy equipartition, because the time scale of the Coulomb interaction is much longer 
than the dynamical time scale. In our model, a; = 10 corresponds to Fr = 2.57. Therefore, 
the threshold x ^ 10 is a border between mild-relativistic and ultra-relativistic shock. If the 
energy of protons can be efficiently transported into electrons only for ultra-relativistic shock, 
the distribution of the break energy can be explained. Hoshino, Shimada (2002) proposed a 
new acceleration mechanism: shock surfing acceleration. This mechanism effectively accelerates 
only electrons in the cases of a strong shock. Such a study might explain the break energy in 
the future. 
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